clear
mu=0.4;
sigma=5;
T=2.1;
lambda=0:0.01:1;
r=length(lambda);

for i=1:r
    y(i)=equas(mu,sigma,T,lambda(i),2,1);       
end

plot(lambda,y)


